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Abstract 

The new model of nn transitions in nuclei based on unitary S'-matrix is considered. 
The \in >-state of nucleus is described by single-particle shell model. The dynamical 
process part is calculated by means of field-theoretical approach with finite time interval. 
The lower limit on the free-space nfi oscillation time r m - m is in the range 10 16 yr > r m - m > 
1.2 • 10 9 s. 
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1 Introduction 



Any information on the occurrence of nn oscillation [1,2] is important in order to discriminate 
among various grand unified theories. The most direct limit on the free-space nn oscillation 
time r min is obtained using free neutrons: r m ; n = 0.86 • 10 8 s [3]. Alternatively, a limit can be 
extracted from the nuclear annihilation lifetime measured in proton-decay type experiments 
(see, for example, Refs. [4-11]). In this case one should calculate the nn transition in nuclei 
followed by annihilation: 

(nucleus) — > (n — nucleus) — >■ M, (1) 
where M are the annihilation mesons. The analogous process in the medium is 

(n — medium) — > (n — medium) — > M. (2) 

The particle oscillations in absorbing matter take place. 

In the standard calculations of ab oscillations in the medium [12-14] the interaction of 
particles a and b with the matter is described by the potentials U a ,b (potential model). \mUb 
is responsible for loss of 6-particle intensity. In particular, this model is used for the processes 

(1) and (2) [4-11]. 

In [10,11] it was shown that one-particle (potential) model mentioned above does not de- 
scribe the processes (1) and (2) and thus total neutron-antineutron transition probability. For 
instance, the total neutron-antineutron transition probability given by the potential model is 
W ~ 1/r (r is the annihilation width of n in the medium), whereas the realistic calculation 
gives W ~ T (see Sect. 5). In the potential model the effect of final state absorption (annihi- 
lation) acts in the opposite (wrong) direction, which tends to the additional suppression of the 
nn transition. So the potential model should be rejected. The S'-matrix should be unitary. 

(For the oscillations in the external field [15,16] the Hamiltonian is hermitian and so the 
absorption is described correctly. The above-given remark holds only for the processes (1) and 

(2) calculated by means of potential model (non- hermitian Hamiltonian). We also note that 
the potential model describes correctly the channel with n in the final state [11].) 

The unitarity of the S'-matrix means that new model should be developed. In [8,9] we have 
proposed the model of the nn transition in medium followed by annihilation which does not 
contain the non-hermitian operators. It is shown in Fig. la. The results are summarized and 
discussed in [17]. In the present paper the nn transitions in finite nuclei followed by annihilation 
(process (1)) are considered. The process model is shown in Fig. lb. The reason is that the limit 
is extracted from the nuclear annihilation lifetime and so one should calculate the process (1) 
and not (2). As we shall see later, the results are the same as for nuclear matter. However, this 
fact is not obvious since the calculations for the processes (1) and (2) are essentially different. 
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A distinguishing feature of the problem under study is the zero momentum transfer in the 
nn transition vertex. Because of this the S'-matrix amplitudes corresponding to the processes 
shown in Figs, la and lb contain infrared divergence. The problem of infrared divergence for 
the particle in the bound state (Fig. lb) is considered for the first time. To gain a better 
understanding of the material, it is desirable to look through the Ref. [17]. 

The paper is organized as follows. In Sect. 2 we formulate the models for the processes 
(1) and (2). In Sect. 3 the diagram lb corresponding to the model with bare propagator is 
calculated. In Sect. 4 it is shown that in the case of S'-matrix problem formulation the process 
amplitude corresponding to the model with bare propagator is singular. The model with the 
dressed propagator is studied in Sect. 5. The results are summarized and discussed in Sect. 6. 

2 Models 

The qualitative process picture is as follows. The free-space nn transition comes from the 
exchange of Higgs bosons with the mass m H > 10 5 GeV [2] and so the subprocess of nn 
conversion is scarcely affected by a medium effects. From the dynamical point of view this is 
a momentary process: r c ~ 1/mjj < 10~ 29 s. The antineutron annihilates in a time r a ~ 1/r. 
We deal with two-step process with the characteristic time T2 ~ r a . 

Thus, the localization of the neutron incide the nucleus does not tend to suppress of the nn 
conversion. This can be also understood using the analogy with the nuclear /3 decay and decay 
of free neutron. It should be emphasized that above-given qualitative process picture does not 
contradict to well-known results on particle oscillations except the absorption channel (see Sect. 
5.2 of Ref. [9]). 

We consider Fig. la. If the antneutron propagator is bare, it contains the infrared singularity 
conditioned by zero momentum transfer in the nn transition vertex. This circumstance changes 
the standard calculation scheme radically. The same is true for the Fig. lb (see Sect. 4). Since 
the process (2) has been considered in details [9,17], we draw analogy with the model used for 
the diagram la. 

We return to Fig. la. The neutron potential U n is included in the neutron wave function 
(unperturbed Hamiltonian) : 

n(x) = VT 1 / 2 exp(— ipx). (3) 
Here p = (e, p) is the neutron 4-momentum; e = p 2 /2m + U n . The interaction Hamiltonian is 

Til = Tinn + 'H, 

T-i - = e -vj/-vl> + H c 

H^t) = J d 3 x-Hi(x). (4) 
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Here l-i n n and H are the Hamiltonians of the nn transition [4] and the n-medium interaction, 
respectively; e nn is a small parameter with e nn — 1/r, where r is the free-space nn oscillation 
time; and are the operators of the neutron and antineutron fields; m n = m n = m. 




Figure 1: nn transition in the medium (a) and nuclei (b) followed by annihilation. 

For the process shown in Fig. lb we take the single-particle shell model of nucleus. The 
initial neutron state is defined by equation of motion: 

(id t - H )n(x) = 0, 

H = -V 2 /2m + U, (5) 

where U is the self-consistent neutron potential. The interaction Hamiltonian is given by (4), 
where T-L is the Hamiltonian of the n-nuclear interaction. 
The neutron state is stationary: 

7ij(t,x) = e~ ie ^nj(x), (6) 

x = (i, x). Here rij(x) and ej are the eigenfunctions and eigenvalues of the Hamiltonian H : 

H nj(x) =e j n J (x). (7) 

The eigenfunctions rij(x) form the complete orthogonal set. 
The Green function of Eq. (5) is defined as 

[i^; - tf (x')]GV, x) = 5 3 (x' - x)tf(f - t). (8) 

Comparing with Fig. la, we see that the neutron plane wave is replaced by the bound 
state wave function (6); the antineutron propagator should be replaced by the Green function 
G(x',x) defined by (8). Both of these processes are described by identical models: The \in >- 
states are the eigenfunctions of unperturbed Hamiltonian. In the case of diagram la, this is 
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the neutron plane wave. In the case of Fig. lb, this is the bound state wave function. The 
interaction Hamiltonian is given by (4). This is the standard formulation of the problem which 
allows to derive the process amplitude directly from interaction Hamiltonian in contrast to the 
model based on diagram technique for direct nuclear reactions [18]. 

In principle, the antineutron propagator can be bare or dressed. In the latter case the calcu- 
lation is standard and simple. The diagram with bare propagator contains infrared divergence. 
Since the corresponding calculations are non-typical, particular attention is given to the model 
with bare propagator. 

We write the general formulas which are used below. Since rij(x) form the complete orthog- 
onal set, the Green function can be represented as [19] 

G(x',x) = -iO{t' - t)^nj(x')n*(x). (9) 

j 

Using the condition of completeness 

]T^(t,x>*(t,x)=5 3 (x'-x), (10) 

3 

one obtains the important relation 




The RHS of Eq. (9) acts as the Feynman propagator transforming the function n,j(x) into 
the point x'\ n,j(x) — > rij(x'). We recall that the basis of plane waves is not used. Equations 
(9)- (11) are valid for any local potential U . 



3 Model with bare propagator 

In this section we calculate the process (1) for the model with bare propagator. The corre- 
sponding calculations are non-trivial: the problem of infrared divergence for the particle in the 
bound state is considered for the first time. First of all we outline a method of calculation of 
the process (2) for the model with bare propagator [8,9,17] because for the process (1) the idea 
of calculation is the same. 

The amplitude corresponding to Fig. la diverges: 

M = e nn G M a , (12) 
1 1 



G = 



en-p 2 j2m-U n + iO 0' 
<fO\T exp(-i J dxU{x)) - 1 1 0n p >= N(2iT) A 5 4 (p f - Pi)M a , (13) 
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Pn = p, e n = e. Here | On p > is the state of the medium containing the n with the 4-momentum 
p — (e, p); N includes the normalization factors of the wave functions. M a is the antineutron 
annihilation amplitude. It contains all the n-medium interactions followed by annihilation 
including antineutron rescattering in the initial state. Due to this the antineutron propagator 
G is bare. 

These are infrared singularities conditioned by zero momentum transfer in the nn transition 
vertex. There is no compensation mechanism by radiative corrections. This is unremovable 
peculiarity. Moreover, for the problem under study the S'-matrix problem formulation (oo, — oo) 
is physically incorrect. For solving the problem the field-theoretical approach with finite time 
interval [8,20] is used. It is infrared free. The problem is formulated on the interval (t/2, —t/2). 
If H = U n =const. (U n is the optical potential of n), the approach with finite time interval 
reproduces all the well-known results on particle oscillations (see Sect. 5.2 of Ref. [9]). 

For the process shown in Fig. lb the zero momentum transfer also takes place and so it 
contains the infrared singularities as well (see Sect. 4). As with Fig. la, we formulate the 
problem on the finite time interval (t/2, —t/2). 

We consider the process (1) on the finite time interval (t/2, —t/2). (The case of S'-matrix 
problem formulation (oo, — oo) is studied in next section.) 

The vector of initial state is 

\0nj>= |0>, (14) 

where | Orij > is the nucleus containing the n in the state j. Since the basis (6) is used, in 
the expressions for the ^-operators ^ n and ^ n the plane waves should be replaced by the 
eigenfunctions rij(x) (Furry representation). Then 

* n \ 0rij>=nj(t,x) |0>, (15) 

where rij(t,yi) is given by (6). We introduce the evolution operator U(t) — 1 + iT(t). In the 
lowest order in H nn the matrix element Tfi(t) is 

f t/2 

T fi (t) = -i< f\ T(exp(-i J ^ dhHjfa)) - 1 \ 0n j >= 

-e n n < f | E T ^) / Me / d 3 x c y n {x c )y n {xc)b+ |0>, (16) 

k=l J -t/2 J 

T k (t) = {-if f t/2 dh... f k " dt k H(t 1 )...H(t k ). (17) 

J -t/2 J -t/2 

In the last multiplier of Eq. (17) we separate out the antineutron field operator \E^(:r fe ): 

H(t k ) = j d 3 x k n(x k ) = j d 3 x k n'(x k )^ n (x k ). (18) 
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Using (15), we obtain 

T fi {t) = ie nn < / I E T fc _i(i) r fc_1 dt k [ d 3 x k H\x k )J k {t k ) I >, (19) 

k=l 7 -'/ 2 7 

J k (h) = f* t2 dt cj d 3 x c < T(* fi (x fc )*fi(a: c )) > n 3 -(x c ). (20) 



fc i ■'-"' 2 

-t/2 

As in the case of plane waves [17], the following relation takes place 



J d 3 x c < T(y n (x k )ty n (x c )) > nj(x c ) = nj(x k ) (21) 



(Schrodinger fields). This relation is analogue of Eq. (11) in the second quantization represen- 
tation. Equation (20) becomes 



J k {t k )=n j {x k ) ( k dt c . (22) 

J -t/2 

As in (15), 

fhj{x k ) |0>= ^ n (x k )d+ |0>= fn(Xk) \ 0nj>, (23) 

where | Ofij > is the n- nucleus containing the n in the state j (with the energy ej). Turning 
back to the Hamiltonian H{t k ) 

J d 3 x k H'(x k )V n (x k ) = H(t k ), (24) 

one obtains 

00 rtk 

T fi (t) = ie nn <f\Y J T k {t) dt c | On, > . (25) 

7 1 J —t/2 



k=l 



Using the formula 

/t/2 rtk-i rtk ft/2 rt/2 ftk-i 

dh... dt k dt c f(h,...,t c )= dt c dh... dt k f(h,...,t c ), (26) 
-t/2 J -t/2 J -t/2 J -t/2 Jt c Jt c 

we change the integration order and pass on to the interval (t, 0). Finally 

T fi (t) = -e nn f dt c <f\ T n {t - t c ) | On, >, (27) 

J 

T n (t-t c ) = Y,Tk(t-t c ) = 
k=i 

00 ft ftk-i 

J2H) h / dh... / dt k H(t 1 )...H(t k ). (28) 

k= \ Jt c Jt c 

Here < / | T n (t — t c ) \ Ofij > is the matrix element of the antineutron annihilation in ft- nucleus 
in a time t — t — t c . 

Equation (27) coincides with (64) of Ref. [9] except that in Eq. (64) of Ref. [9] the matrix 
element < / | T n (t — t c ) \ 0n p > corresponds to annihilation of ft with the 4- momentum p in the 



medium and not n-nucleus. For the problem under study this distinction is inessential. As in 
[8,9], the process (1) probability W(t) is found to be 

Wit) » W f (t) = e 2 nn t\ (29) 

where Wf(t) is the free-space nn transition probability. The result is precisely the same as for 
the nn transition in medium. The lower limit on the free-space nn oscillation time is r^ in = 10 16 
yr. This value is interpreted as the estimation from above. 



4 Infrared divergence 

In this section it is shown that in the case of S'-matrix problem formulation (oo, — oo) the 
amplitude of the model with bare ptopagator is singular: Tfi(t — > oo) ~ 1/0. 

In Eq. (16) we put t — oo. To realize the adiabatic hypothesis, we introduce the multiplier 
exp(— a \ t c |), a > 0. (In the previous section the adiabatic hypothesis has been not used since 
the limiting transition t — > oo was not made.) Then 

ftk 



r k 

Tfi = ~ < f | E T k(°o) / dt c Hnn(t c )e- a ^ | On, >, 
k=i J -°° 

dh... / dt k H(t l )...H(t k ). (30) 



One obtains Eq. (19), where t/2 = oo and 



Since 



J fc (oo) = / d 3 x c f" dt c < T(q n (x k )q n (x c )) > nj (x c )e- a M. (31) 



< T(^ n (x k )^ n (x c )) >= iG n (x k ,x c ) = 

-i0(t k - t c ) ni( x k)n*(x c ) (32) 



(see (9)), Jfc(oo) becomes 

rtk 



J k = [ d 3 x c [ k dt c ^2ni(x k )ni(x c )nj(x c )e aU . (33) 

J J-oo . 

In line with (6), (14) and (27) (see also (44)), n m (t,x) = exp(— i€ m t)n m (x) and 

J k = T,J ^(x^e-^^n^xje^^n^xje-^^e-"^!. (34) 



Taking into account that 



t k 1 pi(ti—£j)tk 

dt c e^- e ^- a ^ = -— -, (35) 



oo I €i — €j — ia 



we get 

1 



Jk = J2~ . _ . _ • ^( X fc) / d 3 x c n*(x c )nj(x c ) 



1 e **** _ 7ij(xife) 

n i( x fe) — — • ( 36 ) 



i 6j — 6j — ia 



As in the case of nn transitions in medium (see Eqs. (53)- (55) of Ref. [9]), the amplitude 
diverges. The value l/(ej — €j) plays the rule of singular propagator. 



5 Model with dressed propagator 

In the model considered above the matrix element < / 1 T n (t — t c ) \ Onj > (see (27)) and ampli- 
tude M a involve all the n-nuclear interactions followed by annihilation including the antineutron 
rescattering in the initial state. In principle, the part of this interaction can be included in the 
antineutron Green function [9,10,17]. Then the antineutron self-energy E is generated. In Eq. 
(36) one should replace 

^ : r 1 ^— ^ I (37) 



€j — €j — ia €j — €j — E — ia 
In this case the amplitude is non-singular and calculation is standard. We consider the process 
(2) for simplicity. As in [10,17], the process probability is found to be 

W d (t) » |f rt (38) 

(E is the parameter), whereas the potential model gives the inverse Independence [4-11] 

W pot (t) » (39) 

The calculations in the framework of unitary model tend to increase the nn transition proba- 
bility: 

w H r 2 

> 1. (40) 



Wpot 4E 2 

The lower limit increases as well. Let r^ in and r pot be the lower limits on the free-space nn 
oscillation time obtained by means of Eqs. (38) and (39), respectively. It is easy to verify that 
[10] 

T min = 7^ T pot- (41) 

If T = 100 MeV and E = 10 MeV then 

4 = 5r pot . (42) 
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For r pot = 2.36 • 10 8 s [7] Eq. (42) gives 




1.2- 10 9 s, 



(43) 



which exceeds the lower limit given by the Grenoble reactor experiment [3] by a factor of 14. 
The parameter E is uncertain. We have put E = ReU n — U n ~ 10 MeV only for estimation (U n 
and U n are the potentials of neutron and antineutron). If E = 0, we come to the model with 
bare propagator. 

If E — > 0, Wd(t) rises quadratically. So r r ^ in is interpreted as the estimation from below 
(conservative limit). 

6 Summary and conclusion 

Result (38) corresponds to the model with non-singular amplitude. Although there is no in- 
frared singularity, this model has essential drawbacks. The model as well as possible suppression 
mechanisms are studied in [9,17]. In present paper the particular attention was given to the 
model with bare propagator since the corresponding calculations are non-typical. 

It is significant that W d (t) rises quadratically as E — > 0. This circumstance should be 
clarified; otherwise the model under study can be rejected. The calculation in the framework 
of the model with bare propagator gives the finite result, which justifies our approach from a 
conceptual point of view and consideration of the model with bare propagator at least as the 
limiting case. In fact this model seems quite realistic in itself [9,17]. In this connection we 
recall the reasons owing to which the approach with finite time interval has been used. 

Since the S-matrix should be unitary, the calculation should be done beyond the potential 
model. However, the S'-matrix amplitude based on hermitian Hamiltonian contains unremov- 
able peculiarity. Moreover, for the problem under study the S'-matrix problem formulation 
(oo, — oo) is physically incorrect [9]. For these reasons the problem is considered on the interval 



On the other hand, if the problem is formulated on the finite time interval, the decay width 
T cannot be introduced since T = J^f^i I S/j(oo, — oo) | 2 /T , T — > oo. This means that the 
standard calculation scheme should be completely revised. 

We point to the important detail of the model which explains the absence of the suppression. 
It is seen from Eqs. (14) and (27) that due to zero momentum transfer both pre- and post-nn 
conversion spatial wave functions of the system coincide: 



Recall that the neutron potential is included in Hq; the n- nuclear interaction H is involved in 
Hi. Since the Hamiltonian H acts on the antineutron, it is turned on following the forming of 



(f,0). 



Onj > sp = \ Onj > 



sp ■ 



(44) 
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the n-nucleus and so Eq. (44) takes place. There is no energy gap U n — U which can lead to 
a very strong suppression of nn transition. The Hamiltonian of the nn transition changes only 
the internal quantum numbers of neutron: H nn \ 0rij>= e nfl \0n,j> [4]. 

Relation (44) explains the result (29). Nevertheless, we view the results of the model with 
bare propagator with certain caution since the process is extremely sensitive to the details of the 
model. Indeed, in the models with bare and dressed propagators the interaction Hamiltonians 
Hj and unperturbed Hamiltonians are the same. The sole physical distinction between models 
is the zero antineutron self-energy in the model with bare propagator. However, it leads to the 
fundamentally different results (see (29) and (38)). This is because the amplitude (12) is in 
the peculiar point. Due to this the problem is extremely sensitive to the value of antineutron 
self-energy E as well as the description of initial neutron state [18] and the value of momentum 
transfered in the nn transitions vertex [17]. 

In conclusion, new model of nn transitions in nuclei based on unitary ^-matrix has been 
considered. Since the results are the same as for nuclear matter, the conclusions are identical 
as well (see Sect. 5 of Ref. [17]): taking into account the result sensitivity to the details of the 
model, the values r^ in = 1.2 • 10 9 s and r^ in = 10 16 yr are interpreted as the estimations from 
below (conservative limit) and from above, respectively. So the realistic limit r m i n can be in 
the range 

10 16 yr > r min > 1.2 • 10 9 s. (45) 

The estimation from below r min > 1.2 • 10 9 s exceeds the restriction given by potential model 
by a factor of 5 and the lower limit given by the Grenoble reactor experiment [3] by a factor of 
14. At the same time the range of uncertainty of r min is too wide. Further investigations are 
desirable. 
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